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Abstract. 

In August 2002, XMM-Newton devoted two full orbits to the observation of 
IE 1207.4-5209 , making this isolated neutron star the most deeply scrutinized 
galactic target of the mission. Thanks to the high throughput of the EPIC in- 
strument, ~ 360,000 photons were collected from the source, allowing for a very 
sensitive study of the temporal and spectral behaviour of this object. The spectral 
data, both time-averaged and phase-resolved, yield one compelling interpretation 
of the observed features: cyclotron absorption from one fundamental (~ 0.7 keV) 
and three harmonics, at ~ 1.4, ~ 2.1 and ~ 2.8 keV. Possible physical conse- 
quences are discussed, also on the basis of the obvious phase variations of the 
features' shapes and depths. We also present deep VLT optical data which we 
have used to search for a counterpart, with negative results down to 1/ ~ 27. 
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1. Introduction 

The X-ray source IE 1207.4-5209 attracted much interest since its early discovery 
(Helfand & Becker 1984) as a bright unresolved source located at the geometrical cen- 
ter of the shell-like radio/X-ray /optical supernova remnant G 296.5-1-10.0 (Roger et al. 
1988). X-ray observations with the Einstein (Helfand & Becker 1984), EXOSAT (Kellet 
et al. 1987), ROSAT (Mereghetti et al. 1996) and ASCA (Vasisht et al. 1997) satellites 
showed a steady flux of ~2xl0~^^ erg cm~^ s~^ (0.3-3 keV) characterized by a thermal 
spectrum. These observations, coupled with the lack of an optical counterpart down to 
V 25 (Bignami et al. 1992; Mereghetti et al. 1996) strongly suggested a neutron star 
nature for IE 1207.4-5209 . This hypothesis was later confirmed by the Chandra de- 
tection of fast X-ray pulsations with period P=0.424 s (Zavlin et al. 2000) and period 
derivative P = {^.OtiD x 10"" s s"^ (Pavlov et al. 2002). 

The Chandra data also unveiled the presence of broad absorption features at energies 
of ^0.7 and ~1.4 keV (Sanwal et al. 2002). The existence of these features was soonafter 
confirmed by an XMM-Newton observation, which showed that the depths and profiles 
of the two lines vary significantly with the rotational phase of the pulsar (Mereghetti et 
al. 2002a, hereafter Paper I). 

Although observationally firmly established, the nature of these lines could not be un- 
ambiguously identified. They were attributed cither to HcII transitions in a "Magnctar"- 
like field B^(1.4-1.7) x 10^"' G (Sanwal et al. 2002) or to atomic transitions in heavier 
elements (e.g. He-like Oxigcn or Neon) in the atmosphere of a neutron star with a more 
conventional magnetic field (Hailcy & Mori 2002). An alternative explanation of the lines 
as cyclotron features was considered to be unlikely by Sanwal et al. (2002). This view was 
criticized by Xu et al. (2003), who interpreted the lines as electron cyclotron resonance 
features originating near the neutron star surface. 

The breakthrough came with a longer XMM-Newton observation which, besides con- 
firming the two phase-dependent absorption lines at 0.7 and 1.4 keV, showed a statis- 
tically significant third line at -^2.1 keV, as well as a hint for a possible fourth feature 
at 2.8 keV (Bignami et al. 2003, hereafter Paper H). The nearly 1:2:3:4 ratio of the line 
centroids, as well as the phase variation, naturally following the pulsar B-field rotation, 
strongly suggest that such lines are due to cyclotron absorption processes in a magnetic 
field of ~8xl0^° G or '~1.6xl0^'' G, respectively in the case of electrons or protons 
features. 
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Thus, among known isolated neutron stars, IE 1207.4-5209 stands out as the only 
one which clearly exhibits X-ray absorption lines. 

Here we present a comprehensive analysis of the data collected during the long XMM- 
Newton observation, already shortly discussed in Paper II, as well as new optical images 
of the field, the deepest available so far, performed with the VLT. 

2. XMM-Newton Data reduction and analysis 

The XMM-Newton observation of IE 1207.4-5209 started on August 4, 2002 and lasted 
two full orbits yielding two uninterrupted time intervals of ~36 hours each. 

The data reported here were obtained with the European Photon Imaging Camera 
(EPIC) instrument, which consists of two MOS CCD detectors (Turner et al. 2001) and 
a pn CCD instrument (Striider et al. 2001), for a total collecting area ^ 2500 cni^ at 
1.5 keV. The mirror system offers an on-axis point spread function of 4-5" FWHM and 
a field of view of 30' diameter. 

While the two MOS cameras were operated in "full frame" mode, the pn camera was 
operated in "small window" mode to allow for the accurate timing of the source photons 
(6 ms resolution). All cameras used the thin filter. The data were processed with the 
latest release of the XMM Newton Science Analysis Software (SAS version 5.4.1). After 
screening the data to remove time intervals with high particle background and correcting 
for the dead time, we obtain a net exposure time of 139.0 ksec for the pn camera and 
197.5 and 197.3 ks for the MOSl and M0S2, respectively. 

The source IE 1207.4-5209 is clearly detected in the three EPIC cameras. In the 
energy range 0.2-5 keV it yields net count rates of 1.466±0.003 cts s~^, 0.380±0.002 cts 
s-i and 0.389±0.002 cts s^^ in the pn, MOSl, and M0S2 respectively. 

The field of IE 1207.4-5209 as seen by EPIC is shown in Fig.[T] Data from the MOSl 
and M0S2 cameras have been merged to produce the image, which has energy-coded 
colors (see caption). 

In the MOS field of view a large 200) number of serendipitous sources are detected. 
Bright emission from parts of the surrounding supernova remnant G296. 5+10.0 is also 
visible. No extended emission is seen around IE 1207.4-5209 or in its immediate vicinity 
(within ^2 arcmin). The radial intensity profile of the target is fully consistent with the 
instrumental point spread function (Ghizzardi 2002). 

To derive the sky coordinates of IE 1207.4-5209 we computed the boresight correction 
to be applied to the default EPIC astrometry. This was done independently for the MOSl 
and M0S2 cameras; the pn data were not used since the small field of view (4x4 arcmin) 
prevented the detection of a suitable number of serendipitous sources. 
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Fig. 1. The XMM-Newton view of the field of IE 1207.4-5209 . Data from the MOSl 
and M0S2 camera have been merged to produce the image. The exposure time is of 
~200 ksec per camera. The target is the bright source close to the center of the field. 
Many (^200) serendipitous sources and the bright emission from parts of the shell of the 
host supernova remnant G296. 5+10.0 are clearly visible. The colors are energy-coded, 
red represents photons in the 0.3-0.6 keV band, green and blue correspond to the 0.6-1.5 
keV and 1.5-8 keV bands respectively 

The positions of the ~200 serendipitous sources detected in the MOS field of view 
were correlated with the Guide Star Catalog II (GSC-II^). After rejecting ambiguous 
matches, the optical brightness and the X-ray spectrum were used as criteria to identify 
the X-ray sources having stellar counterparts. This yielded 6 good reference sources which 
were used to correct the EPIC astrometry. The rms error between the refined X-ray and 
GSC-II positions was found to be of ^1" per coordinate, entirely consistent with the 
expected <1.5" internal EPIC astrometric accuracy (Kirsch 2003). The resulting MOSl 
position of IE 1207.4-5209 is aj20oo = 12/i10to00.91s, 6j2oqo = -52°26'28.8" with an 
overall error radius of 1.5", including the quoted residual uncertainty as well as the 
^ http:/ /www-gsss. stsci.edu/gsc/gsc2/GSC2home. htm 
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absolute intrinsic accuracy (0.35" per coordinate) of the GSC-II. The M0S2 position is 
aj20oo = 12/i10to00.84s, Sj2qqq = — 52°26'27.6" , with an imcertainty of 1.5", similar to 
the MOSl case. The two positions are fully consistent, with a difference well within the 
expected accuracy of the relative astrometry between the MOS cameras. 

In order to obtain an independent measurement on the position of IE 1207.4-5209 , 
we have retrieved from the Chandra Data Archive a public dataset relative to a recent 
(2003/06/15) ACIS Timed Exposure mode observation (20 ksec) of the target. 

Following the Chandra X-ray Center threads to improve the absolute astrometry of 
the standard pipeline-processed archived data^, we used the Aspect Calculator^ to verify 
that the selected observation was not affected by any known aspect offset. 

Then, as in the case of the EPIC data, we used the positions of the serendipitous 
sources in the field to refine the astrometry. We correlated the ACIS and EPIC source 
positions in order to reject any spurious detection in the lower statistic Chandra obser- 
vation. Considering the region within ~ 10 arcmin from the target, we selected 8 secure 
coincidences. Their coordinates were cross-correlated with the GSC-II catalog. Only two 
sources were found to have a match within 3". The boresight correction was found to be of 
^ 0.8", in agreement with the expected absolute aspect accuracy^, with an uncertainty of 
~ 0.35". The best Chandra/ACIS position of IE 1207.4-5209 is a j20oo = 12/il0m00.826s, 
^72000 = — 52°26'28.43" with an uncertainty of 0.6". The more accurate Chandra position 
lies inside the intersection of the MOSl and M0S2 error circles. This gives us confidence 
about the correctness of our analysis and the absence of systematics. 

2. 1 . Timing A n a lysis 

We based our timing analysis on the pn counts extracted from a circular region of 43" 
radius centered on the source position and with energy in the 0.2 - 3.5 keV interval. 
After converting the times of arrival to the Solar System Barycenter, we searched the 
period range from 424.12 to 424.14 ms using both a folding algorithm with 8 phase 
bins and the Rayleigh test. Both methods yielded a highly significant detection at P = 
424. 13076±0. 00002 ms. As in Paper I, the best period value and its uncertainty were 
determined following Leahy (1987) and verified through simulations. 

Comparing the new period measurement of IE 1207.4-5209 with that obtained 
with Chandra in January 2000 (Pavlov et al. 2002) we obtain a period derivative 
f'=(1.4±0.3)xl0~^^ s s~^. This is consistent with the value given in Paper I, but it 
has a smaller error. 

^ http: / / cxc.harvard.edu / ciao /threads/arcscc_corrcction / index.html 
^ http://cxc.harvard.edu/cal/ASPECT/fix_offset/fix_offset.cgi 
" http://asc.harvard.edu /cal/ ASPECT /celmon / 
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Fig. 2. Period history of IE 1207.4-5209 . Circles are the Chandra measurements, crosses 
XMM-Newton 



We note that the P value rests totally on the first Chandra period measurement. 
Using only the 3 most recent values, the period derivative is unconstrained, as is clearly 
seen in Fig. |2 

To study the energy dependence of the pulse profile, we divided the data in four chan- 
nels with approximately 52,000 counts each: 0.2-0.52 keV, 0.52-0.82 keV, 0.82-1.14 keV 
and 1.14-3.5 keV. We verified that an independent analysis of each of these channels with 
the method described above would have resulted in a statistically significant detection 
of the pulsation. The pulse profiles in the different energy ranges (Fig. ^ show a broad, 
nearly sinusoidal shape, with a larger pulsed fraction in the 0.52-0.82 keV band. 

The pulsed fraction in the four energy ranges, defined as the ratio between the num- 
ber of counts above the DC level and the total number of counts, is of (4.0±0.6)%, 
{13.8±0.6)%, (8.3±0.6)% and (9.5±0.6)% (from the lowest to the highest energy range). 

Comparing now the shapes of the light curves, we note that all but the very soft one 
show the peak at the same phase. A phase shift of nearly 90° between the profile in the 
lowest energy range (<0.52 keV) and those at higher energies is apparent. Indeed, a fit 
to the light curves with a sin function yields best fit phases of 0.023±0.026, 0.294±0.006, 
0.245±0.010, and 0.255±0.010 from the lowest to the highest energy range. 
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Fig. 3. Folded light curve of IE 1207.4-5209 in four energy ranges. 



We reanalyzed the data of the December 2001 observation to search for the same 
effect, but the lower statistics hampers a firm conclusion. In fact, using the ~6600 counts 
collected in that observation below 0.5 keV ,the pulsations are only marginally detectable. 
Finally we note that the new observation confirms that no phase shift occurs at ~1 keV, 
as reported by Pavlov et al. (2002) using Chandra data. 

2.2. Spectral Analysis 

To perform the spectral analysis we used events extracted from a 43" radius circle cen- 
tered on the source, selecting PATTERN in the range O-j-4 for the pn and O-j-12 for the 
MOS. As a consistency check, we verified a posteriori that the results of our analysis 
do not change when using only PATTERN (i.e. single pixel) events. Particular care 
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was devoted to the selection of the background regions. In the smah pn field of view we 
excised the source with a circle of 75" radius and we rejected the area possibly contam- 
inated by out-of-time events, or too near to the CCD edges. In the MOS cameras we 
selected a region within the central CCD, excluding both point sources and the diffuse 
emission from the supernova remnant, avoiding the CCD edges, contaminated by Si-K 
internal fluorescence emission. The spectra were rebinned in order to have at least 40 
counts per channel; the pn spectrum was rebinned in order to oversample the instrumen- 
tal energy resolution by a factor of 3. Ad hoc response matrices and ancillary files were 
generated using the SAS tasks rmfgen and arfgen. The spectral analysis was performed 
using XSPEC vll.2 in the energy range 0.3-4 keV. Energies lower than 0.3 keV were 
discarded both in the MOS, owing to calibration uncertainties, and in the pn, owing to a 
strong (possibly electronic) feature in the background at ^0.22 keV. Beyond 4 kcV the 
source is only marginally detected. 

The large number of photons collected in our observation makes statistical errors 
very small. Particular care must be devoted to the systematic uncertainties. The internal 
calibration accuracy of each EPIC detector for on-axis sources is better than 5% (Kirsch 
2003). This yields good quality fits (i.e. <5% residuals at the instrumental edges) for each 
camera (see Fig. 0)), but does not ensure the correctness of absolute flux measurements. 
Indeed, while the cross-calibration between the MOSl and M0S2 cameras agrees within 
5%, differences up to 10% are found between the pn and the MOS, the MOS flux being 
smaller than the pn one below ^1.5 keV and higher above this energy. These effects are 
probably due to remaining uncertainties in the vignetting and CCD quantum efficiency 
functions (Kirsch 2003). 

We therefore adopted the following strategy for the phase-integrated spectra. As a 
first step, the computation of the best fitting model was performed separately for each 
instrument, accounting only for statistical errors. As a second step, to compute reasonable 
confidence intervals on the measured physical parameters of the target, we added an extra 
5% systematic error to each spectral channel. 

A different approach was used for the phase- resolved spectral analysis. For this study 
only the pn data can be used, owing to the MOS slow readout mode. Our aim is the 
description of the relative spectral variations as a function of the pulse phase, which are 
correctly characterized even when the systematic uncertainties are not included. Such 
systematics affect all of the spectra in the same direction (since they are taken with the 
same instrument) and are not a matter of concern when relative differences are studied. 
Therefore, we do not include systematics in the evaluation of the confidence intervals 
for the phase resolved spectral parameters to give a correct description of their relative 
variation, warning the reader that the quoted errors are possibly underestimated when 
the absolute values of the parameters are considered. 
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2.2.1. Phase integrated spectroscopy 

As a first step, we addressed the phase integrated spectrum with respect to our previous 
analysis (Paper II). The improved understanding of the instruments (concerning, e.g., 
the Quantum Efficiency, the Charge Transfer Inefficiency, the redistribution function) 
implemented in the most recent SAS release yields significant differences in the low 
energy portion of the source spectrum, below the O edge (~0.55 keV), with respect to 
our previous analysis. As a consequence, the spectral parameters reported here supersede 
the results of Paper II. 

The spectrum of IE 1207.4-5209 is very complex. Single component (blackbody, power 
law, ...) or double component (blackbody+power law, blackbody+blackbody...) contin- 
uum models alone are totally inadequate. A satisfactory fit requires a model including 
broad absorption features. In the following we will discuss separately the continuum and 
the line components, starting from the pn, which collected ^ 194,000 photons in the 0.3-4 
keV band. 

Single component continuum models fail to reproduce the data (e.g. xt > 5, 96 d.o.f. 
for a simple blackbody, including the line components). The sum of a blackbody and 
a power law yielded a xt ^ 1-6 (94 d.o.f., including the lines), predicting an excess of 
counts at E >2.5 kcV. The best fit continuum curve {xt = 1.15, 94 d.o.f. including the 
lines, see Fig. upper panel) is represented by the sum of two blackbody functions. 
Assuming a distance of 2 kpc (Giacani et al. 2000), the cooler blackbody (hereafter BBl) 
has a temperature kT=0.163±0.003 keV and an emitting radius of 4.6±0.1 km ; the 
hotter (hereafter BB2) one has kT=0.319±0.001 keV and an emitting radius of 800±50 
m. The best value for the interstellar absorbing column is N/f = (1.0±0.1) x 10^^ cm^^. 

We note that a non-thermal, power law spectral component with a luminosity similar 
to that observed for middle-aged pulsars such as PSR B0656-t-14, Geminga and PSR 
B1055-52 (see e.g. Becker & Aschenbach 2002), namely a few 10^*^ erg s"^ in the 0.3-4 
keV range, would not be detectable at the distance of IE 1207.4-5209 . 

The BB parameters differ from those reported in Paper II owing to the updated 
calibrations used here. In particular, the flux at low energy (in the 0.3-0.5 keV range) is 
found to be ~30% higher. 

Four absorption features are clearly seen in the pn spectrum (see Fig 0] and top panel 
of Fig.EJ at the harmonically spaced energies of ~0.7 keV, ~1.4 keV, ~2.1 keV and ~2.8 
keV. 

The different spectral continuum model resulting from the improved calibrations also 
yields a more significant detection of the third and fourth features with respect to Paper 
II. Using a simple gaussian in absorption, we estimated with an F-test that the 2.1 keV 
and the 2.8 keV features have a chance occurrence probability of ~10^^ and ~ 10^"^, 
respectively. 
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In the spectral region encompassing the 2.1 keV absorption feature, two instrumental 
edges, due to Si (1.839 keV) and Au (2.209 keV), are present and, owing to minor 
miscalibrations in the instrumental response, they may give rise to structured residuals 
at a few % level in high-statistic spectra (Kirsch 2003). To address the issue of the 
calibration accuracy in such energy range, we have retrieved from the XMM-Newton 
Science Archive the dataset relative to the observation of a very bright, featureless source, 
the quasar 3C273. We plotted in Fig. ^ (see caption for details on the data analysis) the 
ratio of the data with respect to the best fit continuum model for both 3C273 and 
IE 1207.4-5209 . In the case of 3C273, tiny deviations at the ^6-7% level are seen at 
the expected energies of the edges, with an Equivalent Width of <10 eV. As already 
pointed out in Paper II, the case of IE 1207.4-5209 , is remarkably different. The feature 
at ~2.1 keV represents a 25-30% depletion with respect to the continuum level, while its 
Equivalent Width is ~100 eV. Thus, we rule out the possibility of an instrumental origin 
of the absorption feature near 2.1 keV in the spectrum of IE 1207.4-5209 . 

The large number of photons allows to determine with high accuracy the profiles of 
the two main features. A simple gaussian is inadequate to reproduce the broad dips, 
owing to their asymmetric shape. The combination of two gaussian lines in absorption 
can mimic their profile yielding a statistically good fit {x1 ~ 1.0, 92 d.o.f.). We found 
that a good fit can be also obtained using an asymmetric line profile having the following 
analytic form: 

A{E) = exp[-a(^ - Eo)^ - b{E - E^f - c{E - Eo)^] (1) 

We note that this model does not have a direct physical meaning. Rather, it is used to 
phenomenologically describe the data, since it allows for an easy evaluation of the central 
energy of the features and reproduces their shape with a small number of parameters. 
We therefore adopt this model to fit the two main features at 0.7 and 1.4 keV; the third 
and the fourth features, owing to the lower statistics, can be well reproduced by simple 
gaussian lines in absorption. 

The results of the pn data analysis are reported in Table [2 

We turn now to the MOSl and M0S2 data, encompassing --73,800 and -75,600 
photons in the 0.3-4 keV range. We found that, also in these cases, the best fitting 
continuum model (x^=1.19 for MOSl - see Fig.El middle panel, ^—1A7 for M0S2 due 
to the presence of several small wiggles - see Fig. O lower panel) is represented by the 
sum of two blackbody curves. The temperatures and the emitting regions are consistent 
with the results from the pn camera; the interstellar absorption is found to be somewhat 
higher in the MOS data (N^f = (1.3±0.1) x lO^^cm"^). We ascribe this difference to the 
observed time evolution of the low energy (E<0.5 keV) redistribution function of the MOS 
detectors (Kirsch 2003). This effect is currently under investigation by the calibration 
team; we therefore consider the pn measurement to be more reliable. 
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Fig. 4. Ratio of the data to the continuum best fit models for IE 1207.4-5209 (black 
points) and for the bright quasar 3C 273 (red points). The EPIC/pn observation of 3C 273 
(performed on 2001/06/13; pn operated in Small Window mode with the Medium filter) 
yielded ~62 ksec of good exposure time. PATTERN 0-4 events were extracted from a 
43" radius region centered on the source position. Ad-hoc response files were generated. 
The best fit (xl = 1-29, 1099 d.o.f.) spectral model in the range 0.3-6 keV is a slightly 
absorbed (Nh ^ 1-2 x 10^" cm~^) double power law (photon indices aph i ^ 2.79 and 
aph 2 ^ 1-43; 0.3-6 keV fluxes Fi - 4.47 x 10"" erg cm"^ s^^ and F2 - 8.67 x 10"" 
erg cm"^ s"^). Deviations at the ^ 6 — 7% level are seen in the case of the featureless 
3C273 spectrum owing to residual miscalibrations at the Si (1.839 keV) and Au (2.209 
keV) edges. These artefacts can be reproduced {xl — 1-05, 1093 d.o.f.) by two gaussian 
lines in absorption centered at the edge energies, with equivalent width of ~ 5.5 eV and 
~ 7.7 eV for the Si and the Au edge case, respectively. The ~2.1 keV feature seen in 
the spectrum of IE 1207.4-5209 (see text and Table[21) is markedly stronger, with a >13 
times higher equivalent width, and is totally unrelated to the calibration leftovers which 
are never seen in dim sources. 

Three absorption features are clearly detected in the MOS spectra (see Fig.EJ. The 2.1 
keV feature has a chance occurrence probability of order 10"^ per camera, as estimated 
by means of an F-test. Owing to the lower statistics, the 2.8 keV feature is only marginally 
detected by the MOS cameras and was not included in the model. As in the case of the 
pn, the shape of the two main features was found to be asymmetric and well reproduced 
by the analytic profile described in Equation 1. The parameters of the 0.7 and 1.4 keV 
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data and folded model 




0,5 1 2 

Energy (keV) 



Fig. 5. Fit of the phase- integrated data. The model (double blackbody plus line compo- 
nents) is described in the text. From top to bottom, the panels show data from the pn, 
the MOSl and the M0S2 cameras. In each panel the data are compared to the model 
folded through the instrumental response (upper plot); the lower plot shows the residuals 
in units of sigma. 

features are consistent with the pn results (see Table |5|); however, due to the smaller 
number of photons in the MOS data, some of the parameters decribing the second feature 
are not well constrained. The lower equivalent width of the 2.1 keV feature in the MOS 
should not be a matter of concern since in this region of the spectrum the model for the 
MOS and the pn (including a fourth broad feature at 2.8 keV) are different. 
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Fig. 6. Residuals in units of sigma obtained by comparing the data with the best fit 
thermal continuum model. The presence of four absorption features at ~ 0.7 keV, '--^ 1.4 
keV, ~ 2.1 keV and ^ 2.8 keV in the pn spectrum is evident. The three main features 
are also independently detected by the MOSl and M0S2 cameras. 

2.2.2. Phase-resolved spectroscopy 

As seen in our first EPIC observation of IE 1207.4-5209 , the absorption features are 
phase-dependent (Paper I); in Paper II we showed that the pulse phase variations of 
the spectrum are stronger in correspondence of the features, while the continuum does 
not change significantly. This is evident from Fig. \7\ where the folded light curves have 
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Fig. 7. Folded light curves of IE 1207.4-5209 . The energy ranges have been selected in 
order to put in evidence the impact of the phase variation of the features on the source 
pulsation. The three panels to the left show that in the spectral regions less affected by 
the main features the pulsed fraction is of order '^4-6%. The three panels to the right 
show that the pulsation is much higher (^10%) in the energy ranges where the three 
main features are seen. 



been plotted in six energy ranges where lines are either dominant (right panels) or nearly 
absent (left panels). The pulsed fraction is lower (^4-6%) in the spectral ranges where 
the features are less important, while is definitely higher (~10%) in correspondence to 
the three main features. 

In order to further investigate such an effect, we extracted phase- resolved spectra. 
Following Paper II, we selected the phase intervals corresponding to the peak (phase 
interval 0.40-0.65 with respect to Fig. |21and Cj), the dechning part (phase 0.65-0.90), 
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the minimum (phase 0.90-1.15) and the rising part (phase 0.15-0.40) of the folded Ught 
curve. 

The resulting spectra were fitted allowing both the thermal continuum and the lines 
to vary. For the sake of simplicity, the fourth feature, owing to its lower significance, was 
not included in the model. 

We give the best fit parameters, describing the phase resolved spectra, in Table |21 
Note that systematic uncertainties were not included, as stated in Sect. 12.21 In Fig. |S1 
we plot the residuals with respect to the continuum to show the phase variations of the 
features. 

The results of the phase resolved spectroscopy can be summarized as follows: 

— The two components of the continuum vary slightly (see TableOJ both in temperature 
and in flux with the pulse phase, accounting for ^--^3-5% of the source pulsation. A 
comparison with the values reported in Sect. and an inspection of Fig. [7| clearly 
shows that the phase variation of the features is largely responsible for the observed 
pulsation of the source. 

— the features are strongly dependent on the pulse phase; variations in their width, 
depth (see Table and shape (see Fig. (S)) are evident. The central energy of the 
0.7 keV feature (as derived from the phenomenological model) vary at most of 6%, 
while displacements of the other features are not significant. The 2.8 keV feature is 
marginally detected only during the minimum and the rise intervals. 

— The relative intensity of the first three features varies in phase. 

3. Optical data analysis 

The field of IE 1207-5209 was observed in Service Mode between April and July 2002 
with the 8.2-meter UT-1 Telescope (Antu) of the ESO VLT (Paranal Observatory). 
Observations were performed with the the FOcal Reducer and Spectrograph 1 (FORSl) 
instrument, operated as imager in its standard resolution mode, with a pixel size of 0'.'2 
and a corresponding field of view of 6f8 x 6.'8, large enough to include many reference 
stars for a precise CCD astrometry. Images were acquired through the Bessel V and R 
filters for a total integration time of ^2 and 3 hrs, respectively. The total integration time 
was split in short exposures of 100 s to avoid saturation of two relatively bright stars 
(star A and B - Bignami et al. 1992) close to the position of our target and allow for an 
efficient cosmic ray filtering. Table ^ reports the complete journal of the observations. 
Exposures were acquired under dark time conditions at airmass (secz) '--^ 1.5 and with 
average seeing of ~ 0'.'7 and ^ 0'.'9 in the V and R bands, respectively. Atmospheric 
conditions were good, with little wind and humidity always below 10%. All nights were 
photometric. 
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0.5 1 2 

Energy (keV) 

Fig. 8. Residuals in units of sigma obtained by comparing the data with the best fit con- 
tinuum model (line components removed) for the phase-resolved spectra. The variations 
in width, depth and shape of the absorption features as a function of the pulse phase 
are evident. Structured residuals are observed in the range 0.4-0.7 keV for the "rise" 
spectrum and (more significantly) in the range 0.9-1.2 keV in the "decline", where broad 
emission features could be present. We did not attempt to model these structures, which 
are likely artefacts due to an inadequate description of the main absorption features' 
profile. 
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Standard reduction steps including debiassing and fiatfielding were applied through 
the ESO FORSl data reduction pipeline. Flux calibration was performed using images of 
photometric standards from the Landolt fields (Landolt 1992), yielding extinction^ and 
color-corrected zero-points with an overall accuracy of a few hundredths of magnitude. 
Night-to-night zero-point variations were found to be below 0.04 magnitudes. For each 
filter, images taken in different nights have been registered and averaged using a 3 cr 
clipping algorithm to reject cosmic ray hits. 

In order to register accurately the target position on the FORSl images, we have 
recomputed the image astrometry using as a reference the positions of 20 well-suited, 
(i.e. not extended, not too faint, and not too close to the CCD edges) stars selected from 
the Guide Star Catalogue II. The overall accuracy of our astrometric solution was of 0'.'17 
per coordinate. 

Fig. |5| shows the inner portion of the combined FORSl y-band image centered on 
the target position, with the MOSl, M0S2 and ACIS error circles superimposed. 

A faint object (marked with the two ticks in Fig. O is detected just outside the 
southern edge of the MOSl error circle and showed variability along the time span covered 
by our observations (see caption of Fig.EJ. Its position, in any case, falls more than 2" 
away from the most probable one and we can rule it out as a potential counterpart of 
IE 1207.4-5209 . 

No candidate counterpart is detected in the Chandra error circle (nor in the intersec- 
tion of the MOS ones) down to R^27.1 and V^27.3, which we assume as upper limits 
on the optical flux of IE 1207-5209. 



Date 


Filter 


No. of exp. 


Exposure(s) 


seeing (") 


airmass 


2002 April 22 


R 


43 


4 300 


1.0" 


1.33 


2002 May 10 


R 


4 


400 


0.8" 


1.034 


2002 July 9 


R 


38 


3 800 


0.58" 


1.134 


2002 July 9 


V 


26 


2 600 


0.58" 


1.134 


2002 July 13 


V 


26 


2 600 


0.52" 


1.054 


2002 July 14 


V 


26 


2 600 


0.96" 


1.017 


2002 July 15 


R 


26 


2 600 


0.96" 


1.017 



Table 1. Summary of the optical observations of the field of IE 1207—5209 performed 
with the FORSl instrument at VLT/Antu. 



^ The atmospheric extinction correction was computed according 
to the most recent coefficients measured for Paranal, available at 
|http://www.eso.org/observing/dfo/quality/FORSl/qc/photcoefi^/ 1 
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M0S2 




Fig. 9. Combined F-band FORSl image of the central region of the IE 1207-5209 field. 
The circles indicate the expected source positions computed from the Chandra/ ACIS 
(0.6" radius), MOSl and M0S2 (2" radius) data. The good agreement among the three 
independent positions gives us confidence on the overall correctness of our analysis. No 
sources down to R~27.1 and V~27.3 are seen inside the error circles. The object labelled 
with the two small ticks, visible also in the NTT observation of Bignami et al. (1992), 
falls just outside of the MOSl error circle. Its time-integrated magnitude in the combined 
images is F ~ 25.8 and R ~ 24.5, with uncertainties of ~ 0.1 magnitudes. Its profile is 
consistent with a stellar-like object. Its fiux varied of ~ 0.4 magnitudes in the y-band 
images (from V ~ 26.2 to F ~ 25.7), collected during 3 nights, and of 1.5 magnitudes 
in the J?-band images, collected during 4 nights. It is most probably a background AGN. 

4. Discussion 

4.1. The thermal continuum emission 

The X-ray spectral energy distribution of IE 1207.4-5209 shows a continuum emission 
of purely thermal origin which is well reproduced by the sum of two blackbody curves. 

Thanks to the unprecedented throughput of the EPIC pn camera, we detected with 
high statistical significance the variation of the X-ray continuum with the pulse phase. 
The two blackbody components show a low amplitude (<5%) modulation in both tem- 
perature and flux. This strongly suggests a non-uniform temperature distribution on the 
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neutron star surface. The colder blackbody component (kT^0.16 keV) is emitted from a 
quite large fraction of the surface (R ^ 4.5 km), while the hotter (kT^O.3 keV) is possi- 
bly coming from a heated polar cap (i? ^ 800 m). The use of a more physical contimmm 
model (i.e. a magnetized atmosphere model) could yield lower surface temperatures and 
larger emitting radii, as observed in several cases (e.g. for Vela, Pavlov et al. 2001). The 
phase shift of ^ 90° observed between the peaks of the light curves below and above 0.5 
keV requires an energy-dependent asymmetry in the emission pattern. 

The optical upper limits fall more than two orders of magnitude above the extrap- 
olation of the best fit X-ray blackbody curves. For a distance of 2 kpc, after dered- 
dening for an absorption Ay ^0.6, the upper limit on the optical luminosity Lopt is 
^ 3.4 X 10^^ erg s^^. Using the estimate of the total rotational energy loss obtained from 
the X-ray timing, Ej-ot ~ (7.4±1.5) x lO"^'' erg s^^, we derive an optical emission efficiency 
Lopt/Erot < 4.6 X 10^^. This value is comparable to the optical emission efficiency of 
middle-aged neutron stars like PSR B0656-I-14 and Geminga. 

4.2. The absorption features as cyclotron absorption lines 

The reanalysis of the long EPIC observation, taking advantage of the recent improve- 
ments in the instrument calibration, yielded a more significant detection of the absorption 
features at 2.1 keV and 2.8 keV. The conclusion of Paper II is therefore strenghtened: 
the presence of four absorption features, having central energies very close to the ratio 
1:2:3:4, coupled with their observed phase variation, is naturally explained by cyclotron 
absorption from one fundamental and three harmonics. 

The absorption features are seen during all phase intervals. The absorbing layer is 
therefore surrounding most (all) of the X-ray emitting region. The observed AE/E of 
the features, about O.l-j-0.3 depending on the phase interval, implies that the variation 
of the magnetic field across the line forming region {AE/E AB/B) is very small, con- 
straining both the thickness and the latitudinal extent of the absorbing layer. Assuming 
for simplicity a standard dipole configuration, the radial dependence of the magnetic 
field yields AB/B ~ 3Ar/r. The data therefore constrain the thickness Ar to (0.3-l)rio 
km, depending on the phase interval, where rio is the distance (in units of 10 km) of 
the absorbing layer from the neutron star centre. Moreover, considering the dependence 
of the dipole field on the latitude, B~ (1-1-3 cos^ 6)i, the measured line width implies 
that the absorbing region is limited in latitude to an interval of 15-30 degrees from the 
equator, or 30-50 degrees from the pole, depending on the phase interval. 

4.3. The cyclotron magnetic field and the pulsar slow-down 

The cyclotron interpretation of the absorption features allows for a direct measure of 
the neutron star magnetic field. This can be compared with the magnetic field value 
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estimated from the observed timing parameters, B = (2.6±0.3) x 10^^ G, assuming an 
uniform slow-down due to magneto-dipole braking. 

Assuming that the 0.7 keV feature corresponds to the fundamental cyclotron energy 
for electrons or for protons, the inferred magnetic field would be of 0.6(l+z) x 10^^ G or 
1.2(l+z) X 10^^ G, respectively: z represents the gravitational redshift where the absorp- 
tion occurs. Two hypotheses about the absorbing region position can be explored: (i) in 
the atmosphere, close to the NS surface and (ii) in the magnetosphere, at a few stellar 
radii. 

In the case of electrons, if the cyclotron features are formed close to the surface, the 
fundamental cyclotron energy of ^0.7 keV yields B^8 x 10^" G , assuming a standard 25% 
gravitational redshift. This value is ^ 30 times lower than expected from the observed P 
and P. A "braking problem" arises: some additional torque should be acting in order to 
produce the observed spin-down of the neutron star. A disk from the fallback of supernova 
ejecta could induce such an additional braking if the system is in the propeller regime. 
Debris disks have been invoked (Chatterjee et al. 2000) to account for the properties 
of anomalous X-ray pulsars (AXPs, Mereghetti et al. 2002b). The disk models predict 
an excess in the optical emission (e.g. Perna et al. 2000); indeed, the detection of a few 
AXPs in the IR band (e.g. Israel et al. 2003 and references therein) recently renewed 
attention to the disk hypothesis. Although the actual properties of such systems are 
largely unconstrained, in the case of IE 1207.4 5209 the optical upper limits would 
require any hypothetical disk to be underlmiiinous*' . As a consequence, this scenario 
appears rather unlikely. Alternatively, the cyclotron absorption could take place in the 
magnetosphere, as discussed by Sanwal et al. (2002). Sturner & Dermer (1991) and 
Dermer & Sturner (1994) have suggested that a layer of accreted plasma could form 
in the magnetosphere of a neutron star, at a height of a few stellar radii, supported by 
radiation pressure. An alternative possibility is represented by the presence of a "blanket" 
of j er pairs formed in the closed line region of the neutron star magnetosphere at an 
altitude of ^3 stellar radii (Wang et al. 1998). Although the actual suspended mass, the 
density, as well as the stability of the suspended absorbing layer are highly uncertain, we 
note that in this picture the inferred surface magnetic field would be very close to the 
value expected from the spin parameters. 

In the case of proton absorption, the inferred magnetic field would be at least of 
1.6 X 10^^ G (absorption occurring at the surface). Such a "magnetar" -like field would 
not be free from problems when considering the measured P. The spin-down expected 
from the standard dipole formula would be much larger than observed, 6 x 10~^^ 

^ a rough estimate with a standard, purely dissipative disk model (assuming an inner radius 
equal to the magnetospheric radius, tm ~ 1.6 x 10*^ cm and a rate of material reaching tm 
of ~ 6 X 10^^^ Mq y~^) yields a luminosity ~400 times higher than allowed by the VLT 
observations 
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s s^^. A possibility for solving this problem could be the presence of higher multipole 
components in the magnetic field of the neutron star, as suggested by Sanwal et al. 
(2002). The surface value could exceed lO^** G, while the braking would be driven by the 
dipole component, dominating at large radii. 

4.4. A P problem ? 

The magnetic field inferred from the observed spin-down is hardly conciliable with the 
independent estimates offered by the cyclotron interpretation of the absorption features. 
Furthermore, our refined P value (see Sect. I2.1|l confirms the "age problem" (Pavlov et 
al. 2002): the characteristic age of IE 1207.4-5209 , Tc = (4.7±1.0) x 10^ yrs, is more than 
50 times higher than the age of the associated supernova remnant {tsnr ~ 7 kyrs). 

To reconcile this discrepancy the possibility of a birth period close to the present one 
was proposed (Pavlov et al. 2002; Paper 1). The same born-slow hypothesis was suggested 
to solve similar age inconsistencies in few other cases (e.g. PSR J1811-1925 in Gil. 2-0. 3, 
Kaspi et al. 2001; PSR J0538-)-2817 in S147, Kramer et al. 2003). Such sources represent 
a big problem for massive stars core-collapse theory, since it is difhcult to explain initial 
spin periods as large as a few tens of milliseconds (Heger et al. 2003). 

Of course, all of the above assumes IE 1207.4-5209 to be a smoothly slowing down 
pulsar. If the period evolution of the source is not monotonic, our measurement of the 
period derivative, based on a set of sparse observations, could be wrong. For example, 
glitches could dominate the long-term spin-down of the source, which would appear lower 
than that due to the magnetodipole radiation. The P measurement might also be affected 
by Doppler shift, if the neutron star is in a binary system. In this case, the optical upper 
limits {R > 27.1, V > 27.3), for a distance of 2 kpc and an absorption Ay 0.6, exclude 
any main sequence star as a possible companion to IE 1207.4-5209 , leaving open only 
the possibility of a degenerate object. 

We note that a very recent (June 2003) Chandra observation of IE 1207.4-5209 did 
not help to clarify the issue of the period evolution of the source (see Zavlin et al. 2003). 

4.5. IE 1207.4-5209 : out of the chorus line ? 

An important problem still lacking a solution concerns the uniqueness of the phenomenol- 
ogy of IE 1207.4-5209 in the isolated neutron stars panorama. High quality X-ray ob- 
servations have been performed on several isolated neutron stars of various kind. The 
detection of absorption features has been claimed only in the Soft Gamma Repeater 
1806-20 (Ibrahim et al. 2002,2003), in the Anomalous X-ray Pulsar IRXS 170849-400910 
(Rea et al. 2003) and in the isolated neutron stars RBS 1223 (Haberl et al. 2003) and 
RX J1605. 3-1-3249 (van Kerkwijk 2003). None of these cases is comparable to the phe- 
nomenology of IE 1207.4-5209 . 
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A possible explanation invokes selection effects. IE 1207.4-5209 could be an "uncon- 
ventional", low magnetic field neutron star (< 10^^ G as in the electron, near-surface 
cyclotron scenario). In sources with more typical magnetic fields of order 10^^ G the 
cyclotron energy would lie in the range of the tens of keV, where the surface thermal 
emission is very low and where sensitive X-ray observations have not yet been obtained. 
If the magnetic field of IE 1207.4-5209 , on the contrary, is closer to the typical value of 
10^^ G (as in the magnetospheric cyclotron scenario), the absence of cyclotron features 
in other sources could be explained by a limited lifetime for the absorbing layer. 

An answer to the "uniqueness problem" could come only through a better understand- 
ing of the overall properties of the Central Compact Obiects in supernova remnants, the 
fraternity of IE 1207.4-5209 . These sources are supposed to be the youngest members 
of the radio-quiet neutron stars family (including also Anomalous X-ray Pulsars, Soft 
Gamma Repeaters and Dim Thermal Neutron Stars), but their physics remain elusive. 
We do not understand the lack of radio emission, the lack of X-ray pulsations (IE 1207.4- 
5209 being unique, at the moment, also in this aspect of the phenomenology). The fallback 
of supernova ejecta is possibly playing an important role, driving their multiwavelength 
emission, their spin-down and their evolution, as suggested by Alpar (2001). Deeper 
observations of these sources could shed light to the overall scenario. Only in this per- 
spective we will have an answer to the question whether IE 1207.4-5209 is indeed an 
unique object, or is simply in a standard (transient) phase of the evolution of a young 
neutron star. 
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Table 2. Results of phase-integrated spectroscopy 



pn 



MOSl 



M0S2 



Nh (10^1 cm-2) 
kTsBi (keV) 



R 



sii (km) 
kTsB2 (keV) 



R 



(a) 
SS2 



(km) 



l.OiO.l 1.4±0.1 1.3±0.1 

0.164±0.001 0.165±0.002 0.168±0.002 

4.5±0.1 4.6±0.1 4.6±0.1 

0.319±0.002 0.322±0.002 0.320±0.002 

0.83±0.02 0.83±0.03 0.83±0.03 



El (keV) 

ai 

bi 

Cl 

EWi (eV) 
FWHMi (keV) 
Ea (keV) 
a2 
b2 

C2 

EW2 (eV) 
FWHMa (keV) 

E3 (koV) 
(73 (keV) 
EW:i (eV) 
E4 (keV) 
(74 (keV) 
EW4 (eV) 

^ol-ikev (erg s'l) 
L^'^ (erg s-i) 

xVdof 
dof 



0.68±0.01 



+2 



77; 



-3681^2 



6167 



0.24±0.01 
1.36±0.01 

-3701^° 
488ljf 
66±6 
0.18±0.02 
2.14±0.03 
0.17±0.03 

94±15 
2.83±0.05 



0.11; 



■0.04 



01+20 
»-l-30 



2.24x10"'' 
2.1x10^^ 

1.15 
94 



0.69±0.01 

80±3 
-419±15 
749±40 



1051? 



0.24±0.02 
1.39±0.02 

89l?? 
-36013^° 
2201^°°° 

64±7 
0.19±0.03 
2.12±0.07 



0.13 



+0.10 

'0.04 



48±20 



2.10x10" 



1.19 
154 



0.69±0.01 

82±3 
-412±15 
684±40 

104i« 
0.24±0.02 
1.37±0.02 

96if 
-2201^™ 
50011°°° 

67±7 
0.20±0.03 
2.12±0.08 



in+"-"' 

^■'^'^-0.04, 



45±20 



2.17x10" 



2.2x10^^ 2.2x10^ 



1.47 
151 



" Radius at infinity for an assumed distance of 2 kpc. 
Observed flux. 

Bolometric luminosity for d=2 kpc. 
All the errors are at the 90% c.l. for a single interesting parameter 
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Table 3. Results of phase resolved spectroscopy 

Minimum Rise Peak Decline 



Nh (10^1 cm-2) fixed 
kTssi (keV) 


1.0 

0.157±0.001 


1.0 

0.156±0.001 


1.0 

0.171±0.002 


1.0 

0.164±0.002 


Riti (km) 


4.67±0.03 


4.85±0.06 


4.01±0.03 


4.20±0.04 


kTBS2 (keV) 


0.297±0.002 


0.302±0.002 


0.293±0.002 


0.290±0.002 


(a) /T \ 

RsB2 (km) 


l.OOiO.Ol 


0.98±0.03 


1.04±0.01 


1.12±0.03 


Jil (keV ) 


U.DD4±U.UUd 


U.O/9±0.UUo 


0. <Uo±U.OOO 


U.DD8±U.Uiz 


ai 


y4.4±1.0 


98.2I2.6 


-70 -7+1.9 




bi 

Cl 

EWi (eV) 


-499_5 

O A + 21 

848+20 
1141^ 


-4D4_g 
C\ -r- 1 

762+^4 
1031^ 


-62[J_Q 

424_j7 
921^ 


- 403_9 
7241^1 


FWHMi (keV) 


0.22±0.01 


0.19±0.01 


0.21±0.01 


0.32±0.02 




0.46±0.02 


0.43±0.02 


0.29±0.02 


0.30±0.02 


E2 (keV) 


1.37±0.01 


1.38±0.02 


1.33±0.01 


1.35±0.01 


a2 


147±11 


34.5±0.5 


intf 


63.31°:^ 


b2 

C2 


-690tll 
1087lf^^ 


7 QQ + 0.09 
-I.OO_QQl 

f. 91 q + 0.010 
U.^J.O_Q QQ2 


-981113 
1629lg8® 


-12111 
60.3li;° 


EW2 (cV) 


82±6 


66±8 


54±6 


93±7 


FWHMa (keV) 


0.14±0.02 


0.35±0.04 


0.14±0.02 


0.21±0.02 




0.41±0.02 


0.19±0.03 


0.19±0.03 


0.32±0.02 


Es (keV) 


2.16±0.05 


2.15±0.06 


2.12±0.04 


2.08±0.06 


(73 (keV) 


0.09±0.04 


0.14±0.05 


0.08±0.06 


0.08±0.06 


EW3 (eV) 




1021^5 


28tll 


17111 


(6) 

xVdof 


0.25±0.06 
1.09 


0.29±0.06 
1.20 


0.13±0.06 
1.35 


0.15±0.06 
1.01 


dof 


79 


79 


78 


78 



" Radius at infinity for an assumed distance of 2 kpc. 

Relative line depth = 1. - F {Eii„e) / Fc (E(j„e), where F is the total flux and Fc is the 
of the continuum only 

All the errors are at the 90% confidence level for a single interesting parameter 



